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We consider the quantum measurement properties of a driven cavity with a Kerr-type nonlinearity 
which is used to amplify a dispersively coupled input signal. Focusing on an operating regime which 
is near a bifurcation point, we derive simple asymptotic expressions describing the cavity's noise and 
response. We show that the cavity's backaction and imprecision noise allow for quantum limited 
linear amplification and position detection only if one is able to utilize the sizeable correlations 
between these quantities. This is possible when one amplifies a non-resonant signal, but is not 
possible in QND qubit detection. We also consider the possibility of using the nonlinear cavity's 
backaction for cooling a mechanical mode. 

PACS numbers: 



I. INTRODUCTION 

A number of recent experiments have made use of 
driven microwave transmission-line resonators for sensi- 
tive, near-quantum limited measurements. These include 
measurements of the position of a nanomechanical oscil- 
lator near the standard quantum limit [TJ [2] , as well as 
measurements of single and multi-qubit systems in circuit 
QED setups O |4j . Such experiments use the microwave 
cavity as an "op-amp" type amplifier [5], where the sig- 
nal to be detected (e.g. the position a; of a mechanical 
resonator or the a z operator of a qubit) is dispersively 
coupled to the microwave cavity, meaning that the cav- 
ity frequency depends on the signal. When the cavity is 
driven, the resulting modulation of the cavity frequency 
by the signal leads to a modulation of the phase of the re- 
flected beam from the cavity. By monitoring this phase 
(e.g. via homodyne interferometry) , one has essentially 
amplified the signal. 

Most experiments using a cavity for such dispersive 
measurements and amplification have not exploited non- 
linearities in the cavity- the cavity is just a driven, 
damped harmonic oscillator. The resulting measurement 
and amplification properties of the system are well under- 
stood. In particular, it is known that this system can be 
used for quantum-limited linear amplification, meaning 
that the total added noise of the measurement can be as 
small as allowed by quantum mechanics (see, e.g., Ref. [5] 
for a pedagogical discussion). 

While the linear-cavity regime is certainly useful, it is 
also interesting to consider another possibility afforded 
by microwave circuit cavities: they can be engineered 
to have strong Kerr-type nonlinearities through the use 
of Josephson junctions [BH5]. The resulting nonlinear 
cavity can then be used for amplification in ways not 
possible with a linear cavity. Attention has largely fo- 
cused on using such devices in the "scattering" mode 
of operation, where the signal to be amplified is inci- 
dent on the cavity from a coupled transmission line, and 
where backaction effects are irrelevant. Experiments us- 
ing this mode have realized single-quadrature amplifica- 
tion and squeezing [7J 15141 lj; this operation mode has 



also been the subject of many theoretical treatments, e.g. 
Refs. [BJ [7J HH H2] . For qubit detection, another possibil- 
ity is to use the bifurcation in a nonlinear cavity to give 
a latching- type measurement, where the final dynamical 
state of the cavity depends on the initial state of the 
qubit [131 HH] ; this scheme has also received theoretical 
attention [IB] . 

In this work, we will instead study theoretically the 
quantum measurement properties of a driven nonlinear 
cavity in the operation mode most relevant to experi- 
ments in nanomechanics and quantum information, the 
so-called "op-amp" mode of operation described above. 
Unlike the scattering mode studied in [6j [Jj [12], here 
backaction is indeed relevant, and plays a crucial role 
in enforcing the quantum limit on the added noise: to 
reach the quantum limit, the backaction noise must be 
as small as allowed by quantum mechanics [17] . We 
will focus exclusively on regimes where there is no multi- 
stability in the cavity dynamics (in contrast to the bifur- 
cation amplifier setup). We note that experiments using 
a nonlinear microwave cavity amplifier in the "op-amp" 
mode discussed here have recently been performed. Vi- 
jay et al. have constructed a nonlinear cavity formed by 
a SQUID [18] , and have used this to detect a dispersively 
coupled superconducting qubit [TS] • A recent experiment 
by Ong et al. |20j uses a nonlinear microwave formed from 
a transmission line resonator and a Josephson junction 
to detect a dispersively-coupled qubit; this experiment 
also investigated backaction effects. 

Our analysis focuses on operation points close (but not 
past) the bifurcation in the cavity response, a regime 
which yields extremely large small-signal, low- frequency 
amplification gain. The approach we use is standard: 
we linearize the cavity dynamics about its mean, classi- 
cal value, and use the resulting linear quantum Langevin 
equations to study the noise properties of the cavity de- 
tector. This allows us to asses its ability to reach the 
"op-amp" amplifier quantum limit; a related analysis is 
presented in Ref. [5]. Despite this standard approach, 
we find a number of surprising conclusions that seem not 
to have been appreciated in the existing literature. In 
particular, we show that the nonlinear cavity near the 
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FIG. 1: Schematic figure of a realization of the nonlinear cav- 
ity amplifier. The cavity is formed by an LC circuit contain- 
ing a Josephson junction (energy Ej). The cavity is damped 
and driven by a coupled transmission line; 6j n and 6 ou t denote 
the input and output fields in the transmission. The input sig- 
nal z is a flux which controls the value of Ej, and hence the 
frequency of the cavity ui c . An experimental realization of 
this system is presented in Ref. |18| . 



bifurcation is equivalent to a degenerate parametric am- 
plifier (DPA) driven with a detuned pump [2T] . The value 
of this effective detuning is not an independent param- 
eter, and tends to a universal value as one approaches 
the bifurcation. This mapping allows us to derive sim- 
ple, analytic asymptotic expressions for the cavity's noise 
and gain that are universally valid as one approaches the 
bifurcation point. 

In the low frequency limit, large-gain limit, we find 
that the imprecision noise of the cavity is precisely four 
times what would be expected of an ideal, resonantly- 
pumped degenerate parametric amplifier with equivalent 
gain (c.f. Eqs. ( 45 1 , ([46]) ) . We also show, somewhat sur- 



prisingly, that the cavity's backaction noise at low fre- 
quencies is always given by the same, simple expression 
valid for a linear cavity, an expression which is usually in- 
terpreted as describing the overlap of displaced coherent 
states (c.f. Eqs. (49),([50|). We find that the nonlinear 
cavity amplifier is quantum limited at low frequencies, 
but only if one can make use of the large correlations be- 
tween the backaction and imprecision noises. Such corre- 
lations cannot be utilized simply in QND qubit detection; 
hence, the cavity amplifier misses the quantum limit on 
QND qubit detection by a large factor. 

We also use our approach to study the possibility of us- 
ing the backaction of a nonlinear driven microwave cavity 
to cool a nanomechanical resonator. Near the bifurcation 
point, we find an extremely simple expression for the ef- 
fective temperature of the nonlinear cavity's backaction 



(c.f. Eq. (60 1); surprisingly, the only relevant cavity pa- 



rameter is its damping rate k. We also show that a driven 
nonlinear cavity is far better at cooling a low frequency 
mechanical oscillator than a corresponding driven linear 
cavity with comparable parameters. This latter conclu- 
sion matches what was found by Nation et al. [5], who 
studied backaction cooling by a noninear cavity numer- 



ically over a wide range of cavity parameters, including 
regimes where there is bistability in the cavity dynamics. 
Aspects of cooling and heating using a driven nonlinear 
cavity were also addressed by Dykman [2"2"| . 

The remainder of this paper is organized as follows. In 
Sec. [Tl] we review the basics of how one uses a nonlinear 
cavity as a linear "op-amp" style amplifier, and review 
the formulation and origin of the quantum limit applica- 



ble here. In Sec. Ill we show how near the bifurcation, 



the nonlinear cavity is equivalent to a DPA driven by 
a detuned pump. In Sec. |IV[ we use this mapping to 
derive asymptotic expressions for the cavity's noise and 
amplifier gain near the bifurcation point, and assess its 
ability to reach the quantum limit for small signal fre- 
quencies. Sec. |V| extends this analysis to non-zero signal 
frequencies. Finaly, in Sec. |VI[ we consider the asymmet- 
ric quantum backaction noise of the cavity, and assess 
the possibility of using the cavity for backaction cooling 
a mechanical oscillator. 



II. BASICS OF A NONLINEAR CAVITY 
AMPLIFIER 

A. System Hamiltonian 

The Hamiltonian of a cavity detector with a Kerr-type 
nonlinearity has the general form 

H = H sys + H K = %u) c aSa — HAa'a'aa + H K , (1) 

where uj c is the cavity resonance frequency, A is the Kerr 
constant. We take A > in what follows, as is appro- 
priate for a microwave cavity incorporating Josesphson 
junctions; the results are easily generalized to A < 0. The 
term H K represents the damping (at rate k) and driving 
of the cavity due to its coupling to extra cavity modes 
(e.g. in a microwave circuit, to the transmission line used 
to drive the cavity). Derivations of this Hamiltonian for 
microwave circuits incorporating Josephson junctions are 
presented in many places in the literature, and we do not 
repeat them here (see e.g. Refs. [6"H5I IT2l ITS] ) ; a schematic 
is presented in Fig. [TJ In writing Eq. (JlJ we have assumed 
the relevant case of a high-Q cavity and thus made use of 
the rotating wave approximation to write the nonlinear 
term. We will also be interested throughout in the case of 
a weak nonlinearity, A«k. For clarity, we focus exclu- 
sively on the ideal case where there is no internal cavity 
loss; we also focus on the case of a one-sided cavity. Our 
analysis could be easily generalized to incorporate either 
a two-sided cavity or internal loss (see e.g. Ref. [T2"]). 

Unlike a linear cavity, the nonlinear cavity described 
by Eq. can undergo a bifurcation as a function of its 
parameters from a regime where the average cavity pho- 



ton number 



era) is a single valued function of the 



drive frequency uJd, to a regime where it is multivalued. 
For drive strengths just below the bifurcation threshold, 
n is a single-valued function of drive frequency, but ex- 
hibits a very pronounced slope (see Fig. [2}. This ex- 
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treme sensitivity to cavity frequency makes the cavity an 
extremely sensitive dispersive detector and amplifier in 
this regime. However, it is not a priori obvious whether 
the cavity's noise in this regime is small enough to allow 
quantum limited performance. Answering this question 
is our main goal. 



B. Quantum limit on amplification in the op-amp 
mode of operation 

We focus throughout on the op-amp mode of ampli- 
fier operation, where the input signal to be detected (de- 
scribed by an operator z) is coupled directly to the cavity 
photon number: 



flint = Aa* a ■ z — F ■ 



(2) 



The operator z could represent (for example) the posi- 
tion of a nanomechanical beam (as considered in Ref. [5]) 
or the signal flux applied to a SQUID circuit (as in 
Ref. [E]). As a result of this dispersive coupling, the 
cavity frequency and reflected-beam phase shift become 
z-dependent; one can thus amplify z(t) by monitoring 
this phase. We will focus here on homodyne detection, 
where the output beam is interfered with a classical ref- 
erence beam, and the resulting intensity measured using 
a square-law detector. Letting 6 ut(^) denote the output 
field from the cavity (as defined in standard input-output 
theory [23] [23]), the measured homodyne intensity will be 
described by an operator I(t) 



l = By 



t(*)+ h.c. 



(3) 



where <j> is the phase of the classical reference beam, and 
B is a dimensionless constant proportional to the ampli- 
tude of this beam. Note that / has units corresponding 
to a photon flux. As the value of B plays no role in what 
follows (it is just a scale factor for the output), we set 
B = 1 without loss of generality in what follows. 

We will be interested in weak enough couplings that 
our cavity acts as a linear amplifier. As such, we have a 
linear relation between the input signal and cavity out- 
put: 



(/(f)) = / dt' XlF {t-t'){z{t')), 



(4) 



where xiF{t) oc A is the forward gain of the amplifier, 
and is determined by a standard Kubo formula [5]. 

The amplifier output (i.e. /) will have fluctuations even 
in the absence of any coupling to the detector; these are 
described by the symmetrized spectral density: 



die™* ({I(t), 1(0)} 



(5) 



Again, as we are interested in linear amplification, the 
expectation value is taken in with respect to the state 



of the uncoupled detector (i.e. A = 0). It is useful and 
standard to think of these intrinsic output fluctuations 
in terms of effective signal (i.e. z) fluctuations; we thus 
introduce the imprecision noise spectral density: 



S zz [uj] = S7/[w]/ |xif[w]|' 



(6) 



In the op-amp mode of operation, a second crucial as- 
pect of the amplifier's noise is its backaction. By virtue 
of the the detector-signal coupling in Eq. ([2| , the opera- 
tor F = Aa^a (i.e. the cavity photon number) acts as a 
noisy backaction force on the signal. Extra fluctuations 
in z due to this stochastic force will necessarily increase 
the noise in the output of the amplifier, and are thus part 
of the total added noise of the amplifier. The backaction 
force noise is characterized by a symmetrized noise spec- 
tral density S ff [uj] defined analogously to Eq. 

We thus have that the total amplifier contribution to 
the output noise has contributions from both imprecision 
and backaction noise. It is convenient (and common) to 
think of this total added noise in terms of a noise temper- 
ature T/v [to] : the total amplifier added noise at frequency 
uj is equivalent to the extra equilibrium noise we would 
get by raising the temperature of the signal source by 
Tjv[w] [33] ■ This quantity is relevant no matter what the 
signal, be it the position of a harmonic oscillator or the 
voltage produced by some input circuit; we also stress 
that achieving the quantum limit on the noise temper- 
ature is equivalent to achieving the quantum limit on 
continuous weak displacement detection [5]. 

Minimizing the noise temperature at a given frequency 
requires one to first optimize the signal source's suscep- 
tibility Xzz[u\- This linear-response susceptibility tells 
us how the average value of z changes in response to a 
perturbation which couples to z, i.e.: 



6(z(t)) 



dt' Xzz {t-t'){F(t')). 



(7) 



Optimizing the total added noise over the coupling 
strength and phase of the signal source's susceptibility 
XzzM yields a standard bound on T/v[w] [5]: 



k B T N [u] 

%0J 



> 




S zz [uj}S ff [uj} - [Re(S zF [u])]' 



-ImS zF [ui] 



(8) 



where the inequality becomes an equality for an optimal 
source susceptibility satisfying: 



\XzzM = 

Re XzzH = 

\Xzz[u>]\ y/S zz [w]S FF [uj]' 



S zz [uj]/S ff [uj] 
Re S zF [u] 



(9a) 
(9b) 



We have introduced the correlator S zF which describes 
possible correlations between backaction and imprecision 
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III. BEHAVIOUR NEAR BIFURCATION 



Xi f M 



(10) 



Consider the simple case where the signal frequency 
lo is much smaller than the relevant frequency scales of 
the cavity; we may thus focus on the noise temperature 
in the u> — > limit. Using the fact that there cannot 
be any out-of-phase noise correlations at zero frequency 
(i.e. »S 2 f[0] = Re 5 2 f[0]), the zero-frequency form of the 
fundamental Heisenberg inequality on detector noise [5]: 



S zz [0]Sff[0}- (S zF [0}) 2 >^/4, 



implies that 



(11) 



(12) 



i.e. the added noise amplifier must at least as large as 
the zero-point noise of the signal source [51 [T7] [15] . We 
stress that while the conclusion may appear similar, the 
"op-amp" quantum limit considered here is not identical 
to the quantum limit on the "scattering" mode described 
in the seminal works by Haus and Mullen [25] and Caves 
|26| : the scattering-mode quantum limit does not involve 
backaction. Moreover, an amplifier may reach the quan- 
tum limit in the scattering mode but not in the op-amp 
mode [5]. 

The case where the input signal z is the spin operator 
of a qubit is also interesting. Here, the quantum limit 
on QND qubit detection involves the measurement rate 
Tineas and the measurement-induced backaction dcphas- 
ing rate T v [271 [25]: 

(13) 



r > r 

1 ip _ 1 mcas j 



where for weak coupling: 



T v = 2S F F[0]/h z 
= as = (2§„[0])~ 1 . 



(14a) 
(14b) 



Thus, reaching the quantum limit on QND qubit de- 
tection places more stringent requirements on the de- 
tector than those required to have a quantum-limited 
noise temperature: not only must the quantum noise 



bound of Eq. (11) be satisfied as an equality, but in ad- 
dition, there must be no backaction-imprecision correla- 
tions (e.g. S zF [0] = 0). 

In the following sections, we will calculate the nonlin- 
ear cavity's noise and response functions, and determine 
whether it reaches the quantum limit on its noise tem- 
perature, and on QND detection. We note in passing 
that Ref. [5j also addresses the quantum limit on ampli- 
fication (specifically position detection) using a nonlinear 
cavity in a similar regime to that considered here. Their 
analysis is based on alternative formulation of the quan- 
tum limit which is not equivalent to the one discussed 
here; in particular, they did not address whether the non- 
linear cavity optimizes the quantum noise inequality of 
Eq. (11), or consider its noise temperature as defined in 
Eq. j8jK 



A. Mapping to a degenerate parametric amplifier 

We begin our analysis by using standard input-output 
theory [231 HI] to derive the Heisenberg equation of mo- 
tion for the cavity field, in the absence of any coupling 
to the signal: 



dt 



[a, H s 



(15) 



Here, bi n (t) = bi n e~' lu ' dt + £(i) describes the input field 
incident on the cavity from the transmission line; its av- 
erage value &i n describes the coherent drive applied to 
the cavity at frequency Ud = uj c + A, while £(£) describes 
quantum and classical noise entering the cavity from the 
drive port. Without loss of generality, we take the drive 
amplitude b ln to be real and positive. 

We are interested in driving strengths that result in a 
large average number of quanta h in the cavity, but at the 
same time are not so strong that there is multistability 
in the classical cavity dynamics. It is thus useful to write 
the cavity anhiliation operator a as the sum of a classical 
and quantum part: this takes the form 



a(t) = e 



(16) 



The complex number e l ^ a \/n = (a) is simply deter- 
mined by the classical equations of motion, whereas d 
describes the influence of classical and quantum noise 
(and eventually, the coupling to the input signal). We 
have chosen the phase of the second term in Eq. ( 16 1 to 
simplify the following analysis. From Eq. (15), we find 
that the average cavity photon number \JJi is determined 
by the classical equation 



fi [(k/2) 2 + (2Ah + A) 2 ] = k (6 in ) : 



(17) 



We can now use Eq. ( 15 ) to write an equation for d; 



retaining only leading terms k n > 1 yields a linear 
equation: 



dt 



d,H, 



dpa 



-d-^i(t), (is) 



where 



H, 



dpa 



-hA $d + ih- 



- ddj 



(19) 



Eq. ( 19 ) is simply the Hamiltonian of a degnerate para- 



metric amplifier (DPA) driven by a non-resonant pump, 
where a single pump mode photon can be converted into 
two "signal" mode photons and vice- versa (see, e.g. [24]). 
Here, the classical cavity field a plays the role of the pump 
mode, while the displaced cavity field d plays the role of 
the "signal" mode. The effective parametric interaction 
strength g and effective pump detuning A are given by: 



9 = 2An, 
A = A + 4An. 



(20a) 
(20b) 
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The above mapping of the driven nonlinear cavity to 
a detuned DPA is general, and only relies on S > 1. 
We will be especially interested in operating points near 
the point of bifurcation, as these allow a maximal ampli- 
fier gain. As one approaches the bifurcation the effective 
DPA parameters g, A approach universal values. To see 
this, note first that a standard analysis of the classical 
equations of motions shows that the bifurcation occurs 
at a critical drive amplitude foin : bif satisfying: 



1 



bif 



6^3 A 



(21) 



For 6i n < 6in,bif7 u is a single- valued function of A. For 
frin ^in bifj the slope of n versus A is infinite at a single 
point A = Abifj one finds from Eq. (17 1 



A bif 



V3 

K 

2 

1 K 



2^3 A' 



(22a) 
(22b) 



It thus follows from Eqs. (20) that the parameters of 
the effective DPA attain universal values at the bifurca- 
tion: 



A bi f 



K 
K 

2\/3' 



(23a) 
(23b) 



Note crucially that for cavity operating points near 
the bifurcation, the effective DPA pump detuning A is 
nonzero. As we will see in the next subsection, this will 
have a pronounced impact: the amplified and squeezed 
quadratures of the DPA are not orthogonal. This in turn 
has a significant effect on the noise properties of the non- 
linear cavity detector. 



B. Amplified cavity quadrature 

To appreciate the implications of pump detuning in 
our effective paramp model, we consider the equations 
of motion corresponding to Eq. (19). We will be inter- 



ested throughout in parameter regimes where this effec- 
tive paramp has a photon number gain larger than one; 
this necessarily requires g > |A|. In such regimes, the 
analysis is most conveniently presented by first introduc- 
ing canonically conjugate quadrature operators X and P: 

X = ~{e- ie ' 2 d + e i6 ' 2 (24a) 

P = Z*( e -*«/ 2 d- e "/ 2 #) (24b) 

where for the angle 9 (— 7r/2 < 9 < tt/2) is given by: 

sm6 = A/g. (25) 
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FIG. 2: Average cavity photon number n versus drive de- 
tuning A = u>d — oJc, for various driving strengths b ul . The 
curves show the evolution of the cavity response as one goes 
through the bifurcation: solid curves are for &i n < &i n ,bif , while 
the dashed curve is for &i n > 6i n ,bif. The diverging slope of 
dn/dA near the bifurcation allows for amplification with a 
large gain. 
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FIG. 3: Parametric, zero-frequency photon number gain G = 
G[0] versus drive detuning A, for various driving strengths 
hn < &in. bit (same values as in Fig. The gain is maxi- 
mized near values of the detuning where the slope dn/dA is 
maximal. 



As we will see, the above definition ensures that X is 
the quadrature amplified by the cavity. We also define 
corresponding quadratures X[ n and Pi„ of the operator £ 
associated with noise entering the drive port (e.g. these 
are defined by substituting d —> £ in Eqs. (24)) 



With these definitions, the equations of motion are eas- 
ily solved upon Fourier transforming (see Appendix |A| : 



X[uj] = -y/~K^Xl[u]X in [uj] 

-tan(9( X iH -X2[w]) P„M), (26a) 
P[u>] = -V^XaMPinM, (26b) 



() 



where the susceptibilities X11X2 ar e given by: 



XiM = -iw + k/2 - \Jg 2 - A 2 



X2H 



-iuj + k/2 + \Jg 2 - A 



(27a) 



(27b) 



For the case of a resonant pump (i.e. A = 9 = 0), 
these equations take a simple form and describe the usual 
behaviour of a DPA: as g approaches k/2 from below 
(the parametric threshold), kxi[0] — > 00, kx 2 [0] — > 1, 
and X (P) is the amplified (squeezed) quadrature. By 
considering quadratures of the output field leaving the 
cavity, one finds that the photon number gain for the X 
quadrature is given by: 



G[lo] = \1-k Xi M' 



(28) 



We will refer to G[ui] as the "parametric gain" of our 
system in what follows. For a resonant pump, the am- 
plified and squeezed quadratures are clearly orthogo- 
nal (i.e. canonically conjugate). Note that G[u>] has a 
Lorentzian form, implying that there is only appreciable 
gain for frequency in a bandwidth il^, 



n B = i/xi[o]. 



(29) 



We refer to fig as the "parametric bandwidth" in 
what follows; in the large parametric gain limit, J7s ~ 

The situation is more involved in the case of inter- 
est here, where the effective pump is not resonant (c.f. 



Eqs. (23b|), and hence 9 7^ 0. We still have a parametric 
threshold when g approaches J k 2 /4 + A 2 from below; 
as before, «Xi[0] — >■ 00 in this limit while KX2[0] — > 1. 



It is easy to verify from Eqs. (23) that the parametric 
threshold coincides wi th the cavity bifurcation. It also 



follows from Eqs. (26) that for any pump detuning A, 



X is the amplified quadrature: noise (or signal) incident 
in the X quadrature (i.e. X- m ) only drives the X cavity 
quadrature, and is multiplied by the large susceptibility 
X\ ■ The photon number gain for signals in the X quadra- 
ture continues to be described by Eq. (28 1; as expected, 



this gain diverges as one approaches the bifurcation (see 
Fig. [3]). As one approaches the bifurcation, Eqs. (23) im- 
ply that the angle 9 which defines X takes the universal 
value: 



'bit 



tt/6. 



(30) 



More troublesome when A ^ are the dynamics of the 
P quadrature, the quadrature orthogonal to the ampli- 
fied quadrature. For a non-zero detuning, P is not the 
squeezed quadrature. Noise or signals incident on the 
cavity in the P quadrature (i.e. P ln ) appear both in the 
cavity P quadrature (where it is multiplied by the small 
susceptibility X2), as well as in the cavity X quadrature, 



where it is also "amplified" (i.e. multiplied by the large 
susceptibility xi)- 

To summarize, we have shown that near the bifurca- 
tion, the driven nonlinear cavity of Eq. (JlJ maps onto a 
DPA with a non-zero pump detuning A. This non-zero 
detuning means that the dynamics does not correspond 
to the simple situation of canonically-conjugate ampli- 
fied and squeezed quadratures. As we will see, this lack 
of orthogonality will have pronounced implications on the 
cavity noise properties near the bifurcation. 



C. Coupling to signal and cavity gain 

To complete our mapping of the nonlinear cavity detec- 
tor to a DPA, we need to restore the signal-detector cou- 
pling Hamiltonian and consider the forward gain XiF(t) 
of the system (c.f. Eq.Q). This forward gain tells us how 
strongly the input signal z influences the output homo- 
dyne current, and will not be identical to the parametric 
photon number gain G discussed above. While one could 
calculate xi f M directly using a Kubo formula, it is sim- 
pler here to simply re-derive the equations of motion for 
the cavity field including the coupling to z. Retaining 
only leading-order terms in n, the signal-cavity coupling 
Hamiltonian H lnt in Eq. ^ retains the form H ln t = F ■ z, 
with the generalized force operator F taking the form: 



F = Aa f a~ (VmA) (sin vX + cos vP), (31) 



where 



v= 6»/2 + 3vr/4. 



(32) 



We have dropped a constant term in F which can ab- 
sorbed into the Hamiltonian of the signal source. We see 
that in general, the input signal z couples to both X and 
P, and thus will enter the linearized cavity equations of 
motion as a driving term for both these quadratures. To 
be explicit, one should make the following replacements 
in Eqs. (1261): 



Acosf z[w] (33a) 



X- m [u] -> X- m [oj] 



P ia [cj] -> P m [u] + ( \/— Asiniy ] z[u>}. (33b) 



Note that at the bifurcation, the angle v takes on the 
universal value: 




f bi f = tt/12 + 3tt/4 = 5tt/6 = tt - 



'bif ■ 



(34) 



Having characterized the signal-detector coupling, we 
now turn to the output homodyne current (c.f. Eq. ([3])). 
This current is essentially one quadrature of the cavity 
output field, and may be written: 

= v^fcos^XoutH +shxip h Pout[^]) , (35) 
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where <ph is determined by the phase of the reference 
beam used in the homodyne measurement, and the out- 
put operators are given by the standard input-output 
relations [HEH, e.g. : 



X out [u] = X in [w] + y/nX[w] 



(36) 



It thus follows from Eqs. (26) and rt33| that the linear- 



response gain of the cavity amplifier will have the general 
form: 



X/fH = hAV 2tik(Xi ■ XiM + A 2 ■ X2H) , (37) 



where 



10 r 




Ai = cos iph (cos v + tan 9 sin v) (38a) 
A2 = — smv (cosiph tan# + sinco/j) • (38b) 



IV. AMPLIFIER NOISE IN THE LARGE GAIN, 
LOW FREQUENCY LIMIT 

We are most interested in the properties of our cavity 
amplifier close to the bifurcation, where the parametric 
gain defined in Eq.(28) satisfies G[0] = G ^> 1. In this 
regime, we expect amplification of input signals in a nar- 
row band of frequencies ui < Qb ~ k/VG. We thus 
begin our analysis by considering the amplifier noise to 
leading order in the large parameter G and for uj <C ^s; 
the latter condition allows us to take the zero-frequency 
limit of cavity noise and response functions. We also 
assume the ideal case where the cavity is only driven 
by vacuum noise. Note that it is straightforward to use 
Eq. (17) to determine how G behaves as a function of 



driving strength as one approaches the bifurcation from 
below. Assuming that the drive detuning A is always 
chosen in order to maximize G, one finds that near the 
bifurcation: 



G 



,bif 



,bif 



(39) 



We start with the amplifier's forward gain. In the limit 
G -> 00. Eq. p7| yields: 



X/f[0] 



(40) 



As expected, the forward gain is (to leading order) pro- 
portional to square root of the DPA photon number gain. 

Turning to the cavity output noise, we note that for 
G> 1 and for small frequencies, Eqs. (26) yields that 



the cavity P quadrature is negligible in comparison to 
the X quad rature. As such, we can drop the second term 
in Eq. (35), and treat the homodyne current operator / 
as being proportional to A out oc 
G -> 00: 



X 



G, Thus, for 



I[oj] r~j KCOS(fhX[(jJl]. 



(41) 



FIG. 4: Imprecision noise S Z2 [0], scaled by the imprecision 
noise of an ideal degenerate parametric amplifier, § ZZj ideai ver- 
sus parametric photon number gain G. For large G, the im- 
precision noise is four times the ideal value. Solid blue: in- 
crease G by increasing the drive amplitude bin towards 6i n ,t>if, 
using an optimal A for each drive strength. Dashed red: fix 
(W&in,bif) 2 = 0.995, increase G by tuning A to approach 
the optimal value from below. Long dash - short dash, green: 
same as previous, but increase G by tuning A to approach 
the optimal value from above. 



Looking at Eq. (31) for the backaction force operator 

Thus, to leading 
G: 



F, we see a similar argument holds 
order in G, F is also proportional to A oc 



F[w] ~ y/mAsmvX[w], 



(42) 



Thus, to leading order in G, backaction and imprecision 
noises are perfectly correlated with one another, as they 
only differ by a constant. Their spectral densities will 
simply be proportional to the spectral density of the am- 
plified cavity quadrature, X. 



A. Imprecision Noise 

The leading-order-in-G intrinsic output noise of the 
amplifier (i.e. noise in the homodyne current at A = 0) 
thus follows easily from Eq. (41 ) and (26a) (see Appendix 
\An. In the low frequency limit, we have: 



Su[0] 



cos 2 ip h ■ G(l + tan 2 6). 



(43) 



The two terms in the last factor represent two dis- 
tinct physical contributions to the output noise. The 9- 
independent term arises from A-quadrature input noise 
(Ai„) being amplified and appearing in A. In contrast, 
the term proportional to tan 2 9 is a direct consequence of 
the non-zero effective pump detuning A. The resulting 
non-orthogonality of amplified and squeezed quadratures 
causes P-quadrature input noise (Pi n ) to also be ampli- 
fied and appear in the amplifier output I. We thus see 
that A / causes the output noise to be larger than 



what would be expected for a resonantly-pump DPA with 
equivalent photon number gain G. 

Combing the above expression with Eq. ( |40| ) for the 
gain, we find that the low- frequency imprecision noise 
(c.f. Eq.© ) near the bifurcation (i.e. G — > oo) is given 
by: 



§ zz [0] ~ S zzMeal [0]-h imp (e,v), (44) 



where 



$zz, ideal 



f h\ \ 

\AnA 2 J 

1 + tan 2 



(cos v + sin v tan ( 



(45a) 
(45b) 



Here, S ZZj idcai is the imprecision noise of an "ideal" DPA 
in the large gain limit. By "ideal" , we mean a DPA which 
was pumped on resonance, and where the input signal z 
only drives the amplified X quadrature (i.e. the angle 
v in Eq.(31| would be zero). hj mp (9, v) > 1 describes 
the increase of S zz due to the fact that our nonlinear 
cavity does not realize a DPA in this ideal fashion. The 
numerator of hi mp describes the extra output noise due 
to the non-resonant effective pump, as discussed after 
Eq. ( 43 1 . The denominator describes the reduction in 



gain coming from the fact that the signal drives both the 
X and P quadratures. 

Finally, we can further simplify our result by using the 
fact that near the point of bifurcation, the effective DPA 
parameters approach universal values (c.f. Eq. (23)). To 
leading order, we can simply replace 9 by its value at 
the bifurcation #bif = tt/6. We thus obtain our final 
expression for the imprecision near the bifurcation: 

S zz [0] ~ 4S , ZZj jdeal- 



(46) 



In large gain limit, the imprecision of the nonlinear cavity 
amplifier is a factor of four times what would be expected 
from a theoretically ideal degenerate parametric amplifier. 
The behaviour of the imprecision noise relative to the 
ideal value is shown in Fig. |4j 



Backaction noise and backaction-imprecision 
product 



From Eq. (42 1, we see that to leading-order- in-G, the 
low-frequency backaction noise spectral density 5ff[0] 
will just be proportional to the low- frequency output 
noise spectral density 57/ [0]. Using the universality of 
the DPA parameters near the bifurcation, we find that 
in the G — > oo limit: 



Sff[0] 



lA 2 n 



-G. 



(47) 



We see the backaction diverges as the parametric photon 
number gain G; this is a simple consequence of the fact 




100 
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FIG. 5: Backaction-imprecision product S zz [0]5Vf[0] (scaled 
by h 2 G) versus parametric photon number gain G, demon- 
strating the universal scaling predicted for large G. Individual 
curves correspond to the same parameters as in Fig. [4] 
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FIG. 6: Backaction-imprecision correlations, as measured 
by SzzF [0] / y/Szz [0] Sff [0] versus parametric photon number 
gain G. As expected, the curves all tend to the universal 
value of 1 as G — > oo. Individual curves correspond to the 
same parameters as in Fig. [4] 



that our dispersive coupling unavoidably leads the signal 
z to be coupled to the amplified cavity quadrature X. 
The full expression (valid for arbitrary G) is not too un- 
wieldy and is given in the Appendix as Eq. (A14). Note 



that similar spectral densities for a nonlinear cavity were 
calculated using a linearized Fokker-Plank approach in 
Ref. |16j and (in the classical, high temperature regime) 
Ref. [HI- 

Combining our results, we see that near the bifurca- 
tion, the backaction - imprecision product will be much 
larger than the minimum value of Ti 2 /4 allowed by quan- 
tum mechanics. In the large G limit, we have: 



Sff[0]S zz [0] ~^h 2 . 



(48) 



Figure [5] shows the scaling of S,f_f[0]S' zz [0] versus para- 
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metric gain G as one approaches the bifurcation by either 
tuning the drive detuning A or the drive strength 6; n ; the 



universal asymptotic behaviour described by Eq. (48 1 is 
clear. 

The above result implies that if one cannot make use 
of backaction-imprecision noise correlations, one is very 
far from having a quantum limited device. In particular, 
near the bifurcation the nonlinear cavity detector cannot 
be used for QND qubit detection: in such an experiment, 
the backaction dephasing rate will be a factor G ^s> 1 
larger than the minimum rate dictated by quantum me- 
chanics (c.f. Eq. (13l). Note that the situation is very 
different for a linear cavity: there, as long as one drives 
the cavity on resonance, the 5*f_f[0]5 2Z [0] product attains 
the minimum possible value of h~ /4 [5J [3U] . 

It is tempting to think that by simply changing the 
cavity operating point slightly, one could achieve a situa- 
tion where the input signal is only coupled to the cavity 
quadrature P, and thus avoid the problematic diverging 
backaction found above. From Eqs. (31 ) and (32), we see 



that this would require an operating point for which the 
angle = it/2. However, from Eqs. (25) and (27), this 



in turn implies that the cavity would have no paramet- 
ric gain: G — 1. Thus, one cannot solve the problem of 
large backaction by simply changing the drive detuning 
without simultaneously getting rid of the amplifier gain. 



Comparison with linear-cavity backaction 
formula 



For a linear cavity, one can directly connect the back- 
action noise spectral density at zero frequency to how 
strongly the average cavity amplitude (a) changes in re- 
sponse to a change in the signal. One finds: 



d(a) 

1a 



(49) 



This elegant result was first derived in Ref. [3T] in the 
case where z is a spin operator for a qubit; S*f_f[0] in this 
case is directly proportional to the qubit dephasing rate 
(c.f. Eq.(14a)). Heuristically, it expresses the fact that 



the backaction disturbance of the measurement is directly 
related to the distinguishability of cavity states associ- 
ated with different values of the input signal. A small 
change in the input signal causes a small displacement of 
the coherent state describing the cavity. Eq. ( 49 ) implies 
that the backaction dephasing T v (and hence Sff) is di- 
rectly determined by the overlap between this displaced 
coherent state and the original coherent state describing 
the cavity. 



One would not expect Eq. ( 49 ) to apply in general 



to our nonlinear cavity detector, as now the intracavity 
state corresponding to a given fixed value of the input 
signal is not a coherent state, or even a pure state [32] . 
This is a direct result of the squeezing and amplification 
of the cavity noise that occurs as one approaches the bi- 
furcation. However, if one is far from the bifurcation, 



these effects should be minimal, and one might expect 
Eq. ( 49 1 to remain valid. This idea was recently put for- 
ward in Ref. |20j . and derived within an approximation 
which neglects noise squeezing of the cavity. Our ap- 
proach fully accounts for the squeezing of the intracavity 
fluctuations, and allows us to test the general validity of 
Eq. ( 49 1 . Surprisingly, we find that this expression ex- 



actly captures the full backaction noise, even close to the 
bifurcation: 



SW,lin[0] = SW[0]. 



(50) 



Here 5ff[0] is the full expression for the backaction 
noise spectral density that follows from Eqs. ( [26] ), (see 
Eq. | A14[ )). We see that despite the fact that the cavity 
is not in a coherent state or even a pure state, Eq. (49) 



remains valid for the nonlinear cavity amplifier; that this 
should be so is by no means a priori obvious. 

Ref. [20 also suggests that the nonlinear cavity de- 
tector reaches the quantum limit on QND detection 



(c.f. Eq. (13)), implying that the backaction-imprecision 
product SffSzz attains its minimum possible value of 
h 2 /4. In contrast, we find that the backaction noise Sff 
(and hence backaction dephasing rate) is factor G ^> 1 
larger than the quantum limited value (c.f. Eq.(48) and 
Fig. [5]) . The discrepancy here arises from the fact that 
Ref. piO] does not explicitly calculate the measurement 
rate (i.e. 1/S ZZ ) for a specific, optimized cavity readout 
scheme, but rather assumes that it also be given (up to a 
prefactor) by overlap expression in Eq. (49). This would 
imply the measurement imprecision noise S^fO] scales 
like 1/G in the large-G limit. In contrast, we explicitly 
consider homodyne detection of the cavity output. We 
find that the imprecision noise (and hence measurement 
rate) are independent of G in the large gain limit, in 
agreement with Ref. [18]. This is a simple consequence 
of the fact that the nonlinear cavity's parametric gain 
amplifies both the signal and the vacuum fluctuations 
driving it by the same factor of yG. 



D. Quantum limit on the amplifier added noise 

While in the low- frequency, large gain limit, the nonlin- 
ear cavity system cannot function as a quantum-limited 
QND qubit detector, it may nonetheless be a quantum 
limited linear amplifier (i.e. have the minimum noise tem- 
perature In allowed by quantum mechanics). This dif- 
ference stems from the fact that when used as an am- 
plifier (in the "op- amp" mode), one can take advantage 
of correlations between backaction and imprecision noise 
by tuning the susceptibility of the signal source (e.g. in 
a voltage amplifier, the source impedance). 

To leading order in G and at low frequencies, we have 
shown that the backaction and output noise operators are 
proportional to one another, implying perfect correlation: 



S zF [0] ~ J S zz [0}S FF [0}. 



(51) 
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Fig.[6]shows the behaviour of these correlations versus G, 
where G is tuned is various ways; the asymptotic, perfect 
correlation behaviour is clear. 

Turning to Eq. ([8| for the optimized noise temperature 
Tjy, we see that perfectly correlated backaction and im- 
precision noises do not contribute. This implies that our 
leading-order-in-G analysis is insufficient to determine 
whether Tjv is quantum limited: this analysis only tells 
us that there is no order- \[G term in T/v- To determine 
whether the quantum limit is reached near the bifurca- 
tion, one must go beyond leading order expressions, even 
though we are interested in the low frequency limit. Such 
an analysis is straightforward though tedious; details are 
presented in Appendix|Aj Obtaining the cavity noise cor- 
relators and forward gain exactly from Eqs. (26) with 



no large-G assumption, we find that at zero frequency, 
the nonlinear cavity detector always optimizes the quan- 
tum noise inequality of Eq. (_?_?| ) (i.e. it is satisfied as an 
equality). As such, the minimal low- frequency noise tem- 
perature given by Eq. Q is indeed the quantum limited 
value of Tiuj/2. We stress that this result is completely 
independent of the choice of homodyne phase iph. 



E. Utility of backaction - imprecision correlations 



see that XzzM is purely real if one is far from resonance, 
i.e. \oj — uim\ 3> 7. Thus, the nonlinear cavity detector 
is ideally suited to applications where one is interested 
in non-resonant position detection. For example, stan- 
dard interferometric gravitational wave detectors require 
sensitive position detection of a test mass in the free- 
mass limit, i.e. lom — > |33) . In this case, as long as 
lo 7, one always has a non-resonant situation, and 
Xzz is real. For such frequencies, the nonlinear cavity 
amplifier would be able to achieve a quantum-limited 
noise temperature. In contrast, if one used a detector 
with S Z F — in this regime, the noise temperature is 
at best a factor (Jm/i 3> 1 larger than the quantum 
limited value. The utility of using correlations between 
backaction and imprecision noise is well-known in the 
gravitational wave community [34] , though it is not usu- 
ally discussed in terms of the general noise temperature 
language used here. 

Finally, we note that if the input signal z was a volt- 
age, and we think of our cavity amplifier as a voltage am- 
plifier, the requirement that the input susceptibility be 
purely real to optimize the noise temperature translates 
into requiring a signal source with a purely imaginary 
source impedance [S]. 



As always, achieving a quantum-limited noise tempera- 
ture is not simply a question of having an amplifier which 
saturates the fundamental quantum noise inequality of 
Eq. ( 11 )- one also needs to optimally tune the susceptibil- 
ity x^zM of the signal source (i.e. the source impedance). 
This optimization results in two conditions, c.f. Eqs. Q. 
The magnitude condition (c.f. Eq. (9a)) can always be 
achieved by an appropriate tuning of the signal - detec- 
tor coupling A] it corresponds to properly balancing the 
relative contributions of backaction and imprecision noise 
to the total added noise. In contrast, the phase condition 



(c.f. Eq. (9b)) cannot be achieved by simply tuning A. 
It corresponds to optimizing \zz [H to optimally make 
use of in-phase backaction - imprecision correlations de- 
scribed by Re S z p. 

Consider the nonlinear cavity detector in the low- 
frequency, large gain regime considered above. We f ound 
that it has a maximal value of correlations S z f, Eq. (51 ). 



Eq. ( 9b ) then implies that reaching the quantum limit on 
the noise temperature requires Im XzzM — 0. This is in 
sharp contrast to the more common situation where S z f 
vanishes, and the optimal source susceptibility Xzz must 
be purely imaginary. 

This has interesting consequences. For concreteness, 
consider the case where our input system is a mechanical 



oscillator and z 
given by: 

Xz 



represents a position, XzzM is simply 



— 1/m 



lUJ"f 



(52) 



Here, com is the resonance frequency of the mechanical 
oscillator, m is its mass, and 7 is its damping rate. We 



V. AMPLIFIER NOISE AT NON-ZERO 
FREQUENCIES 

It is straightforward to extend our analysis to describe 
the amplification of signals with frequencies co that are 
non-zero, but still small enough that the parametric gain 



G[(J\ > 1. It follows from Eq. (28J) that in the G[0] > 1 
limit, this requires co < n B ~ n/y/G\0\. Simple ana- 
lytic expressions are easily obtained in the limit where 
G[0] — s- 00 while co/^Ib stays finite. To leading order 
in G[0], one finds (as expected) that the photon number 
gain G[uj] and forward gain x/fM have a Lorentzian 
frequency dependence on a scale set by Qb- Letting 
co = lo/Qbi we have: 



G[co] = 



G[0] 



1 + w 2 



r 1 a 
XIFM = y 



Xif[0] 



(53a) 
(53b) 



In the same limit, we find that the imprecision noise is 
frequency independent, whereas the remaining correla- 
tors also decay with frequency on a scale set by fig: 



S zz [ui] 
§ff[u] 

S zF [i0] 



Szz[0] 

Sff[0} 

1+LO 2 

s zF [o] 

1 + iCo 



(54a) 
(54b) 

(54c) 
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It immediately follows that for finite frequencies, the 
noise temperature behaves as: 



k B T N [u} 



We thus see that at finite frequencies, the reduced noise 
temperature 2feflTjv /hu rapidly increases as a function of 
frequency from the quantum- limited value of 1; in partic- 
ular, it is already much greater than one for frequencies 
small enough to not appreciably reduce the gain. The 
leading correction at finite lo comes from the imaginary 
part of the noise cross-correlator S z f [w] ■ As discussed ex- 
tensively in j5] , such out-of-phase backaction- imprecision 
correlations cannot be taken advantage of by simply tun- 
ing the susceptibility of the source; as a result, their ex- 
istence represents unused information, and thus leads to 
a departure from the quantum limit. In principle, such 
correlations can be utilized via feedback techniques. 



VI. BACK-ACTION COOLING 

We have seen in the preceding analysis that near the 
bifurcation point, the backaction noise of the nonlinear 
cavity amplifier diverges; this prevents quantum-limited 
amplification unless one can make use of noise correla- 
tions. In this section, we change focus somewhat and 
consider the specific case where the input signal z is the 
position of a mechanical resonator. In this case, the large 
backaction of the nonlinear cavity may actually be use- 
ful: it has the potential to strongly cool the mechanical 
resonator towards its quantum ground state. 

The topic of backaction cooling has received consider- 
able attention in the optomechanics and electromechanics 
communities [35] . It has been shown that the backaction 
of a linear cavity dispersively coupled to a mechanical 
resonator can be used to ground-state cool the mechani- 
cal resonator if one is in the so-called good cavity limit, 
where the mechanical frequency u>m is much larger than 
the cavity damping rate k |361 I37j . This regime has been 
exploited in recent experiments with linear microwave 
cavities 138 , 32] and optical cavities [40H42] . 

In the opposite regime of a low-frequency mechanical 
resonator (u>m *C k), cooling using a linear cavity is still 
possible, but the lowest achievable temperature is on the 
order of Tea ~ n/ks- A crucial parameter is the back- 
action damping (or optical damping) rate 7ba : this is 
the enhanced damping of the mechanical resonator re- 
sulting from a net energy loss to the driven cavity. The 
cooling power of the cavity backaction will be directly 
proportional to 7ba- I n the low frequency regime, a sim- 
ple classical linear-response argument yields that for a 
mechanical resonator dispersively coupled to a cavity, 
7b a oc dn/dA [37] 03]. As discussed in Ref. [8], one 
thus expects that a nonlinear cavity will be capable of 
much stronger backaction damping than a linear cavity, 



given the enhanced slope of the cavity response curve (c.f. 
Fig.§. 

A large backaction damping is not, however, in itself 
enough to ensure good cooling. One needs that the cav- 
ity acts as a source of cold damping for the mechanical 
resonator. Thus, one must also consider the effective 
temperature of the backaction noise, Tba- Ref. [8] exam- 
ined this quantity numerically; in contrast, our approach 
allows us to obtain simple analytic expressions in the in- 
teresting regime where one is near the bifurcation point 
and the parametric gain G[0) ^> 1. 

Our analysis is based on the unsymmetrized back ac- 
tion noise spectral density, defined as 



'ffM 



dt(F(t)F(0)). 



(56) 



The symmetrized noise considered in previous sections 
is given by Sff[w] = (Sff[w] + Sff[— k>])/2- The fre- 
quency asymmetry of Sff[w] describes the asymmetry 
between emission and absorption of energy by the cav- 
ity; a standard perturbative calculation shows that it di- 
rectly determines the backaction damping of the mechan- 
ical resonator [5]: 



7baM 



Irahio, 



-(Sff[+u] - Sff[-u]). 



(57) 



where m is the oscillator mass. Using the linearized- 
Langevin approach described in previous sections, we find 
a particular simple form for 7ba near the bifurcation, in 
the limit where um/^b remains constant as the para- 
metric gain G — > oo: 



1 A 2 n G[0] I Afn ^ . f . 

7BA ~ v1^ 2 i + (*) ! = 7i h^ G[u)l (58) 

In the u>m limit, this reproduces the classical expres- 
sion 7b a oc dn/dA, while for non-zero frequency, we see 
that the backaction damping decays rapidly on the scale 
of the parametric amplification bandwidth The full 
expression (valid even for small G) is given in Appendix 
[A] It is instructive to compare this result for 7ba against 
the corresponding expression for a linear cavity, in the 
relevant limit wm <C and for an optimized detuning 
[3"B1 [3"7] . As expected, one finds that the nonlinear cav- 
ity's 7ba is enhanced by a factor of the parametric gain 

As already discussed, we must also consider the effec- 
tive temperature TbaM of the backaction, a quantity 
which is in general frequency dependent and is defined 
as [5] 



1 A 2 - 



exp 



kid 



fcsr BA [w] 



S 



FF-W 



s 



FFl 



(59) 



Using our linearized Langevin approach, we find a partic- 
ularly simple asymptotic expression for the Bose-Einstein 
factor riBA associated with TbaM i n the large-G limit 
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relevant near the bifurcation: 



Acknowledgements 



1 + 2«ba [w] = coth 



Tilo 



2fc B r B AH 



1 + 3(w/k) 5 
v3(w/k;) 



(60) 



The full expression for «ba is given in Appendix [Xj The 
expression for tibaM is remarkably similar to the corre- 
sponding expression for a linear cavity [351 [37] ■ In partic- 
ular, the relevant frequency scale is k, and not the much 
smaller scale set by the parametric bandwidth One 
thus finds that in the low-frequency limit ujj -C k: 



T BA [0] 



2^3' 



(61) 



In contrast, in the low- frequency limit, the effective back- 
action temperature of an optimally-driven linear cavity 
is k/2. 

Thus, the effective backaction temperature of our non- 
linear cavity near the bifurcation only differs by a nu- 
merical prefactor from that of a linear cavity. For low- 
frequencies, the final oscillator temperature T osc is given 
by [37]: 



IqTq + Tba^ba 

70 + 7BA 



(62) 



Here, 70 is the oscillator damping resulting from its in- 
trinsic (i.e. non-backaction) sources of dissipation, and 
To is the temperature of this bath. We have thus es- 
tablished that for a low frequency mechanical resonator, 
the nonlinear cavity is a far better way to cool than the 
linear cavity. One has a much greater backaction damp- 
ing rate, as well as a slightly smaller effective backaction 
temperature. 



VII. CONCLUSIONS 

In this paper, we have given a theoretical treatment of 
the quantum measurement properties of a driven nonlin- 
ear cavity used as a linear detector or amplifier. By using 
the equivalence between this system near its bifurcation 
point and a degenerate parametric amplifier driven by a 
detuned pump, we were able to give a relatively simple 
description of the physics. We find that quantum limited 
amplification is indeed possible, but only if one is able 
to make use of the large correlations between backac- 
tion and imprecision noises. Such correlations are ideally 
suited to position detection of a mechanical system far 
from resonance; however, they cannot be utilized in QND 
qubit detection, and hence one is far from reaching the 
relevant quantum limit on this task. We also examined 
the possibility of backaction cooling using this system, 
demonstrating that the nonlinearity is particularly use- 
ful in the case where one wants to cool a mechanical 
resonator whose frequency « k. 
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Appendix A 



1. Mapping to the Detuned DPA 



Using Eqs. (18 ) and (19), the equation of motion for 
the displaced cavity annihilation operator d takes the 
form: 



d = (-«/2 + iA) d + gd) - %/«£(*)• 
Introducing the canonical quadratures 



1 

71 



1 1 

—i i 



(Al) 



(A2) 



and defining x- ln ,pi n to be the corresponding quadratures 
of the noise operator £, the equations of motion take the 
form: 



d_ fx 
dt \P 



= M 



Here, M is the matrix defined as 

-A 



M = 



9-2 
A -(9 



(A3) 



(A4) 



Eq. ( A3 ) can be conveniently solved by first diagonal- 



izing M. The only subtlety is that due to the nonzero 
effective drive detuning A, M is non-Hermitian; as a re- 
sult, its eigenvectors are not orthogonal to one another. 
Defining as per Eq. (25), we let: 



cos(0/2) sin(0/2)' 
sin(0/2) cos(0/2) 



(A5) 



denote the matrix whose columns are the eigenvectors of 
M. One then has: 



M = V 



(Xi[0])" 






(X2[0])" 



V 



(A6) 



where the eigenvalues of M are just the inverses of the 
susceptibilities Xi[0],X2[0] defined in Eq. (27). 

The rotation to the quadratures X and P introduced 
in Eq. ( 24a ) can now be written as 



cos (0/2) sin (0/2)' 
-sin (0/2) cos (0/2) 



(A7) 



The form of M makes it clear that X as defined in 



Eq. (24a) is indeed the amplified "eigenquadrature" of 
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the cavity: it corresponds to the first eigenvector and 
eigenvalue of M. In co ntrast, the orthogonal quadrature 
P defined in Eq. \2Ah\ does not correspond to an eigen- 
vector of M. 

Finally, Fourier transforming the equations of motion 
Eq. (A7) using the convention: 



We have used the value of the angle v given in Eq. ( 32 1 



A[oj] ee / dtA(t)e~ lult . 

J — oo 

and making the above rotation, we find: 



(iwl + TMT -1 ; 



X[u] 
P\u] 



X- m [uj] 

Pin[u] 



(A8) 



(A9) 



where the identity matrix is ljj = 6ij. Solving for X[uS\ 
and P[oj] directly yields Eqs. (26). 



2. Back Action Force 

Using the definition of the backaction force operator F 
given in Eq. (31 ) and solutions to the cavity equations of 
motion, Eqs. (26), we find: 



F[u] = F x [oj] X in [oj] + F p [oj] 4M, 



(A10) 



where 



F x [lo] ee — A\J 2nn sin v\\ M , 

-Fp[w] ee — A\/2nK[(x2[^] — XiM) sin v tan 9 

+ X 2 Mcosz/], (Alia) 



and the angle v is defined in Eq. ( 32 ) 



The unsymmetrized force noise spectral density 



Sp p [to] defined in Eq. ( 56 ) can be written in terms of 
F[u] as: 

2tt5(u + u}')SffM = (^[w]f[w']). (A12) 



We can thus use Eq. (A10) to calculate SffM if we 



know the correlation functions of the input noise opera- 
tors X- la , P m . From standard input-output theory, and 
our assumption that £(i) describes vacuum noise, one 
easily finds: 

(X in [uj]X in [io']) = tt5(lo + lu') (Al3a) 

(AnMAnK]) = v6(u+u/) (Al3b) 

(X in [uj]P in [uj']) = inS{uj + uj'). (A13c) 

Explicitly computing the symmetrized spectral density 
Sff[w] = (Sff[u] + Sff[— w])/2, we find: 



S FF [w] = 4A 2 nn x 



(A14) 



k 2 + 6g 2 + Auj 2 - 2g 2 (cos 29 + 4 sin 6) 
, (k 2 + 4w 2 ) 2 - 8g 2 cos 2 9{k 2 - 4w 2 - 2g 2 cos 2 9) 



We stress that this expression only involves our initial lin- 
earization of the dynamics, and does involve any further 
assumption of being close to the bifurcation. In the limit 
where one approaches the bifurcation (i.e. KXi[0] — > 00, 
7r/6), one obtains the asymptotic form given in 



Eq. (47) 



One can use Eq. (A14) to verify that 5ff[0] is in- 
deed related to the derivative of (a) with respect to A 
as per Eqs. (49) and (50). This is easily done using 
(6) = v / ^ e ^°: where n is given by Eq. ( |17| , and the 
phase <p a is given by: 



tan < 



2A + 4An 



(A15) 



(as follows from the classical equations of motion). 

Finally, we note that the above results are easily gen- 
eralized to fi nite t emperature. The input noise correla- 
tors in Eqs. (A13) and <§ffM are simply multiplied by 
(1 + nth), where n t h is a Bose-Einstein factor evaluated 
at the temperature of the incident thermal noise, and the 
frequency of the cavity. 



3. Imprecision Noise 

The intensity I[u>] of the homodyne measurement is 
given by Eq. ( 35 1 . Using the solutions to the equation of 
motion in Eq. ( 26 ) we obtain 



= I x [u] X- m + I p [u] P„ 



(A16) 



where 



V / KCOS( f fh(l — K XiM) (Al7a) 
v / K(sin(^/ l (l - KX2M) (A17b) 
-KCOs^tanfl^M - XiM))- 



It is now a straightforward exercise to compute the spec- 
tral density Sjj[uj] from Eq. (A16) and Eqs. (A13), in 
complete analogy to our calculation of S FF [w[- One 
finds that this output noise completely symmetric in fre- 
quency: Sii[u} = Sii[—uj], and thus Su[u>] = Su[u]. 
The imprecision noise spectral density <S zz [w] then fol- 
lows using Eqs. ^ and (37). 



4. Imprecision - backaction correlation 

The symmetrized imprecision-backaction noise corre- 
lator SjfM may be written: 



S/fH = - {Si F [uj] + Sif[-w]*) ■ 
where the unsymmetrized correlator is given by: 
2n5(uj + uj')S IF [uj} = (I[uj}F[uj'}} . 



(A18) 



(A19) 
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We may thus calculate Sif[uj] using Eq. (AlO), (A16) 
and (A13); dividing by xif[u] as given in Eq. (37) then 



yields the desired correlator S z f[lu]. 



5. Cooling 



Using Eqs. (All ), on finds that the full expression for 



the asymmetric- in- frequency part of Sff[w] is given by: 



Sff[+u] — Sff[— w] 



(A20) 



64A 2 ngn (sm9 - 1) 



(k 2 + 4lu 2 ) 2 - 8g 2 cos 2 (k 2 - 4w 2 - 2g 2 cos 2 0) ' 

This expression then directly gives the backaction damp- 
ing via Eq. (57 1. 



Combining the above expression with Eq. (A14) for 
Sff[w], we can find a general expression for u-qa M > the 
effective temperature of the backaction expressed as a 
number of quanta (c.f. Eq.(60l). We have: 



2n B A[w] 



6g 2 + 4w 2 - 2g 2 (cos 20 + 4 sin t 



8gu(sm0- 1) 



(A21) 



We stress that Eqs. (A20l and (A21) do not involve an 



assumption of being near the bifurcation. 



We derived Xi.f[ w ] in the main text by accounting for the 
coupling to z in the cavity equations of motion, resulting 
in Eq. p7|. 



In general, an amplifier may also have a reverse gain 
Xf/[w]; this describes how signals coupled to the output 
operator I could affect the average value of the backac- 
tion force operator F [5]. In general, reverse gain is unde- 
sirable, as it implies that measuring the detector output 
(by coupling to it) can lead to enhanced backaction fluc- 
tuations. The forms of the fundamental quantum noise 



inequality of Eq. (11) are also modified in the presence 
of reverse gain. 

To show that the reverse gain of our cavity amplifier 
vanishes, we make use of the equation [5]: 



Xif[u] - Xf/H* = -{i/fi)[Si F [u} 



Sif[-c 



(A23) 



Using the solution of the cavity equations of motion 
to calculate SjfM, and using the expression for x/f[w] 
from Eq. ( 37 1 , Eq. ( A23 1 directly yields that there is no 



reverse gain at any frequency 



6. Reverse Gain 



The forward gain in our system was defined in Eq. Q , 
which upon Fourier transforming takes the form: 



Xfi[u] = 0. 



(A24) 



= Xif[u]z[w] 



(A22) 
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